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In a previous experiment imbibition fronts mov- 
ing up in a paper sheet, that drifts down into a water 
tank at speed v, were found to have an average height 
width w ~ hP K , and a temporal correla- 
tion function C{t) = v~ K f(tv (6t+K - )/f) ). We simulated the 
phenomenon using a pipe network model and reproduced 
the rich scaling behaviors with all exponents f2, n, 6 t and 
(3 in agreement with the experimental values || . Another 
study based on a phase-field model leads to predictions 
on exponents that are unfortunately incompatible with 
the experimental ones Q . 

The front propagation exemplifies two related dynam- 
ics, namely mean flow and roughening. Assuming a 
unique time scale controlling both dynamics, an expo- 
nent identity /3 = 0.(6 1 + /c)/(f2 + 1) was derived. In- 
serting experimental values into the right-hand side, we 
get [3 = 0.52 compared with the measured value 0.56. 
Using numerical estimates, f3 = 0.54 is obtained to be 
compared to the numerically measured value 0.63. The 
identity is also supported by a constancy in the relative 
degree of roughening for simulated fronts in stationary 
paper sheets 0|. 
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FIG. 1. Log-log plot of characteristic time t\ and ti re- 
spectively for flow and roughening against interface height h 
from experiment (open symbols, right and bottom axes) and 
simulation (solid symbols, left and top axes). 

The unique time scale assumption and hence the iden- 
tity are the main objections by Dube et al. We now 
give further evidences by directly extracting relevant time 
constants from experimental and numerical data in Refs. 
The mean flow is characterized by t\ = h/v. For 


roughening, we observe that C (t) can be nicely fitted by 
V2w[l + (t/t 2 y al3 ]~ 1/a at a = 3.5 and the relaxation 
time t 2 hence follows. Results are plotted in Fig. 1. The 
unique time scale assumption means that t\ and t 2 are 
different representations of a unique time scale for the 
front propagation dynamics. It requires t\jt 2 to be a 
constant of order, but not necessarily exactly unity. In 
contrast, if flow and roughening were distinct dynamics 
||, t\ and t 2 could easily differ by a factor of, say, 10 
or 100. Using the experimental results, t\jt 2 w 0.82 on 
average for various values of h. It is close to unity and 
is a strong experimental support of our view. A further 
numerical support is that ti/t 2 ss 2.0 roughly from our 
simulation. The very similar magnitudes of t\ and t 2 
both in experiment and simulation independently, in ad- 
dition to their approximate proportionality should not 
be causally dismissed as accidental. Instead, its impli- 
cation of a unique time scale can be a helpful hint for a 
complete theory of the experimental observations. 

The discrepancy between the ratios 0.82 and 2.0 de- 
serves some caution, though it can be due to finite size 
effects, which affects prefactors more severely than expo- 
nents. In Fig. 1, variations in the slopes of the fitted 
lines for either the experimental or numerical values re- 
flect slight deviations from t\ <xt 2 . These are alternative 
representations of deviations of the estimated exponents 
from our identity discussed above. They have been at- 
tributed to finite size effects which are particularly sig- 
nificant in the simulations 
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